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We developed a method of precise isotope labeling to visualize the continuous growth of graphene by chem¬ 
ical vapor deposition (CVD). This method allows us to see in real time the growth of graphene monocrystals 
at a resolution of a few seconds. This technique is used to extract the anisotropic growth rates, the formation 
of dendrites, and the dependence on adsorption area of methane on copper. We obtain a physical picture of the 
growth dynamics of graphene and its dependence on various parameters. Finally, our method is relevant to other 
CVD grown materials. 


INTRODUCTION 

Visualizing the growth of materials has been of interest for 
a long time. This can be as simple as watching water crys¬ 
tals grow on a windshield or more complex such as observing 
the growth of proteins in a liquid environment. In most cases 
this is done in-situ via optical means. However, at the atomic 
level this becomes more difficult and other techniques need 
to be used such as LEED or scanning probe techniques IH . 
Typically, crystals grow in an “unfriendly” environment, with 
very high temperatures or extreme pressures, which doesn’t 
permit an in-situ observation. In our work we developed a 
method to continuously image and recreate the crystal growth 
ex-situ, i.e., after the crystal is fully grown. We have applied 
this technique to graphene crystals, which is a fascinating sys¬ 
tem by itself, but the method is generalizable to many other 
materials. This is achieved by varying the relative ratio of the 
isotopes continuously during the growth process. The isotope 
concentration will therefore provide a unique label determin¬ 
ing when the material was grown. This method is similar to 
carbon dating, where the isotope concentration of radio active 
carbon decays over time and therefore allows the dating of or¬ 
ganic materials O. The difference here, are the time scales 
involved (seconds) versus hundreds of years and the spatial 
dependence, which allows us to correlate the position with 
time. 

The ability to visualize the growth of atomic materials has 
important implications, particularly with the recent develop¬ 
ment and potential of hybrid materials beyond graphene Gllll, 
where the atomic layering and growth history is an important 
factor. Eor graphene, in particular, this method not only gives 
new insights on the growth dynamics, like the formation of 
dendrites, of nucleation sites or additional layers as discussed 
here, but the varying isotope method can also be used to la¬ 
bel certain areas of the crystal without changing its electronic 
properties. Other applications include the possibility to tune 
the nuclear spin density Q, the ability to do phonon engineer¬ 
ing reduce thermal conductivity due to phonon disorder 
Q, determine chemical and vibrational properties of different 
layers, |[8l, increase sensitivity to determine specific phonon 
modes (91, and others (ToHEI. 

Here we provide new results organized along three main 
themes: continuous labeling and determination of the sen¬ 


sitivity of the method, which can provide close to hundred 
independent labeling or time steps; dynamic growth or the 
making of movies and time snapshots of the growth, in par¬ 
ticular the fractal growth of graphene; growth rates and the 
dependence of the growth rates on anisotropy, dendricity and 
growth diffusion area. 

Graphene is a two dimensional material which has attracted 
a great deal of attention for its unique physical and elec¬ 
tronic properties ifTsMHl . Originally isolated by mechani¬ 
cal exfoliation of graphiteO, more scalable methods neces¬ 
sary for production of large scale graphene layers and indus¬ 
trial applications include chemical vapour deposition (161 [TtI 
and graphitization of SiC lfT^ . Chemical vapour deposition 
(CVD) of graphene sheets on commercial polycrystalline cop¬ 
per foils is a leading candidate for the production of large 
high quality graphene films. The main challenge in produc¬ 
ing high quality graphene by CVD is controlling the nucle¬ 
ation density of graphene crystals which in turn limits the 
size of individual crystals. Grain boundaries where differ¬ 
ent misoriented graphene crystals come together are known 
to reduce the mechanical strength and electronic properties of 
graphene (1911201. As a result it is crucial to be able to produce 
large monocrystal graphene sheets. Since the pioneering work 
of Li et al. in graphene CVD on copper foils a great deal 
of work has been done in determining the optimal conditions 
to produce large single crystal graphene and understanding the 
role of various growth parameters in the CVD process (211 . 
This has led to large gains in the size and quality of single 
crystal graphene fiakes which can be produced (2^ [23l up to 
mm and cm scales. Analyzing and understanding the growth 
dynamics of CVD graphene, is one of the main contributions 
of this work, where we start by describing a new labeling 
method. 


HIGH RESOLUTION CONTINUOUS LABELING METHOD 

In general, the possibility to label a material locally is of 
importance in a number of applications. Indeed, it allows to 
encode information, without modifying the basic properties. 
In the case of graphene, this can be particularly useful in areas 
where transparency, organic composition or low mass is im¬ 
portant, such in labeling electronic devices or biomolecules. 
The labeling is done by using different relative concentrations 
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FIG. 1. a) Raman map showing 2D position for the single sample containing different ratios of graphene. Scale bar is 10 /xm b) 

Representative Raman spectra for the different isotopic ratios showing the 2D peak, c) We observe a linear shift in position dependent on the 
amount of vs. present during growth. The position of the 2D peaks shown as a function of fraction. Solid grey line indicates a 
linear curve fit to the experimental data, with the line thickness representing a 95% confidence interval. The dashed line indicates the theoretical 
curve given by equation 1. 


of isotopes, which can then be used for encoding or trans¬ 
parent labeling. The detection is achieved by Raman spec¬ 
troscopy, which is used to characterize graphene m and 
other graphitic materials ||25]| . Here we use a combination of 
precisely controlled isotopic methane flow and and 
Raman mapping to probe the growth dynamics of graphene 
crystals in much greater detail than is possible when using 
only a binary on/off isotopic variation method ll26ll . 

Several samples were prepared using pure ^^CH 4 , pure 
^^CH 4 and various mixes. In the carbon isotope graphene 
we observe a downshift in Raman peak position by a factor of 
approximately ^12/13 as a result of the change in the mass 
of the carbon atoms ||9l[T7l. A sample was prepared where the 
ratio of methane was modified stepwise throughout 

the growth in 10% intervals to give ratios of 9:1, 8:2, etc. A 
Raman map of this sample is shown in figure where we 
observe bands of graphene with different 2D peak positions. 
Plotting the position of the 2D peak for each of these bands 
we find and excellent agreement with the predicted linear be¬ 
haviour where the position of the peak is a function of the 
average mass of carbon atoms in a given region. Taking x to 
be the fraction of and Acc = we 

have 

^2D{t) = ^2 dC^C) + XAUJ ( 1 ) 

The fit is done to a linear function with ^ 2 dC~‘^C) and 


^ 2 dC'^C) as free parameters and gives values of 2686.7 cm“^ 
and 2580.6 cm“^ respectively. 

The theoretical line plotted in figure [T]: shows the expected 
result taking uj 2 dC'‘^C) as 2687 cm“^. The excellent agree¬ 
ment quantifies the precision with which the isotopic ratio of 
our graphene during CVD growth and extract this value by 
Raman spectroscopy of our samples. The average standard 
deviation for the peak position for a given concentration is 
cr = 1.4 cm“^. This translates into a labeling resolution of 

Ri =cr/Acj- 1.3*10“^ (2) 

The high resolution allows the labeling of 1/Ri 80 regions 
independently. Moreover, most physical properties are not af¬ 
fected by the labeling, except for Raman process. This can be 
used for a variety of applications, such as encryption, variable 
nuclear spin concentration, or even isotope concentration de¬ 
tection. In what follows, we will use it in order to gain insight 
into the growth mechanism of graphene. 

DYNAMIC IMAGING 

In situ observation of the CVD growth process is techni¬ 
cally demanding Ill and as a result the growth mechanisms of 
the chemical vapour deposition process are not fully under¬ 
stood. Most of the understanding of the CVD growth process 
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FIG. 2. Visualizing the growth of a graphene crystal: Snapshots of the graphene growth showing the combined intensity of the 2D peak as 
a function of time. Intensities are scaled to differentiate between layers 1 and 2. Last square: Raman maps of a dendritic graphene crystal 
indicating the Raman shift of the 2D peak and associated time from the peak position to time correspondence. Only the layer 1 peak position 
is shown. 


is the result of ex situ analysis of samples prepared by CVD 
across a range of conditions l[2T1l . The use of Isotopic methane 
is a useful tool that has been used to characterize the growth 
of polycrystalline sheets (TTl [26l and graphene single crys¬ 
tals. Typically this has been done by using alternating flow 
of and methane to produce graphene which displays 
a tree-ring like growth pattern ifTTll . Here, several graphlocon 
samples were prepared with the ratio of to varying 
linearly as a function of time. The sample shown in flgure|^ 
was first annealed in hydrogen at 1073°C for 4 hours to pre¬ 
pare the surface and then graphene was deposited during a 
400 second growth phase with a combined methane flow of 
1.2 seem. Since the Raman peak position is proportional to 
the isotopic concentration (figure [^, we can then extract the 
expected 2D peak position at any given time using eq. 1. 

The growth conditions employed for this sample were de¬ 
signed to give an incomplete growth with dendritic graphlo- 
cons Eli. A Raman map was taken with grid spacing of 900 
nm. A spectrum was taken at each point on the grid with a 


collection time of Is. The position,width and intensity of the 
2D peak were then extracted. We observe that the 2D posi¬ 
tion is lowest at the center of the graphlocon and increases as 
we move towards the edges. This is consistent with a growth 
that begins with primarily methane, with the concen¬ 
tration increasing with time. The analysis could equally be 
performed using the G or other Raman graphene peak, how¬ 
ever the 2D peak was chosen for the combination of large am¬ 
plitude and high Raman shift, which enable the position of 
the peak to be extracted with greater precision and lower col¬ 
lection time relative to other peaks. We also observe in the 
center of the graphene crystal a region of bilayer graphene, 
characterized by a double Raman peak structure. The growth 
times for each graphene layer are calculated from the position 
of the two independent 2D peaks observed. In fact, we ob¬ 
served two types of bilayer growth, either Bernal stacked or 
not, sometimes even in the same crystal. In the case of Bernal 
stacking, the 2D Raman peaks of each layer hybridize 1^ to 
form one peak corresponding to the average mass between the 
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two layers, while in the non-Bernal stacking, the two Raman 
peaks identify the frequencies corresponding to the respective 
atomic mass of each layer, with a weak enhancement due to 
graphene enhanced Raman scattering |[29i . This is indicated 
in figure where the bilayer is shown in orange. In the last 
square of figure we show only the contribution of the first 
Raman peak corresponding to the first layer. 

This linear shift in peak position allows us to extract an 
absolute time scale for the progression of the crystal forma¬ 
tion. We can visualize the growth of the graphlocon as a func¬ 
tion of time by considering the points in the Raman map with 
U 02 D < ^ 2 D{t). In figurej^we show snapshots of the graphlo¬ 
con at 8 different points in time. In the supplementary material 
the full movie of the graphene growth is shown. 

GROWTH RATES: ANISOTROPY, DENDRICITY, AND 
DIFFUSION AREA 

Its clear from the isotope labelling that the graphene crystal 
grows from the center outwards. This is consistent with previ¬ 
ous studies on CVD graphene growth |[26||30l[3Tl. In some of 
these studies the crystal symmetry was also determined and 
various shapes have been observed, including four-fold and 
six-fold symmetries, depending on the growth conditions and 
substrate (5 ^ . However, the local growth rates have not yet 
been analyzed. Using our data from the method described in 
section 3 we are able to extract the growth rates and directions 
at every point in the crystal. 

Anisotropy 

The growth anisotropy is not affected by the growth stage 
and is constant over time during the growth and independent 
of the size of the crystal. The radial growth is indeed linear 
with radius (see figure 3c) and therefore time, since we chose 
X ^ t. Great care was taken to calibrate the and 
methane fiows, so that the total methane fiow is constant over 
time to allow for a precise determination of the time depen¬ 
dence. To analyze the anisotropy in more detail, we compare 
the growth rate in different directions for a weakly dendritic 
graphene crystal, shown in figure [^. This sample was pre¬ 
pared with a combined Ci74 fiow rate of 1.2 seem and a H 2 
fiow rate of 80 seem, during a 3 minute growth phase. The 
highest growth rates are clustered along 6 arms spaced by 
roughly 60°. By taking profiles along the fastest and slowest 
directions we observe that the average growth rate is roughly 
doubled along the fast profiles when compared to the slow 
ones, 9.0 /im/min vs 4.7 /im/min. 

The growth rate was extracted for growths performed 
across a range of methane gas fiows. Comparing the cal¬ 
culated partial pressure of methane gas with the growth rate 
we observe a clear linear relationship. Figure shows the 
relationship. In all cases this data was measured along the 
fastest growth direction, and includes samples prepared us¬ 


ing the technique discussed in this article as well as 

pure samples where the growth rate is calculated from the 
diameter of the graphlocon and the growth time. To probe 
whether this relationship holds within the same single crys¬ 
tal, we synthesized sample F using a variable overall methane 
fiow as well as shifting the ratio and the same range 

of growth rates as those extracted from individual growths was 
obtained. 


Fractal shape and dendridty 

In ref. it was shown that fractal graphene crystals can 
be synthesized using CVD. Fractal shapes can lead to an en¬ 
hanced broadband efficiency in antenna designs lf33l and the 
corresponding increased edge to surface ratio, can lead to an 
enhanced reactivity for chemical processes at the edge 1341 . 
Pushing the anisotropy further could even lead to the synthe¬ 
sis of quasi one dimensional arms or graphene nano-ribbons. 

A large dendritic sample was synthesized using low pres¬ 
sure and low gas fiow rates. The copper foil was sonicated 
in acetic acid for 30 minutes prior to growth in order to sup¬ 
press graphene nucleation and allow larger domains to form. 
A high resolution Raman map shows details of the formation 
of the dendritic arms of the graphene crystal. The growth rates 
were extracted along the primary arm of the crystal as well 
as several dendrites and shows the dendrites forming concur¬ 
rently with the main arm. The average linear growth velocity 
is extracted by fitting the distance along each profile to the 
time corresponding to the 2D peak position at each point on 
the profile. We see that the main arm has a velocity of 18 
yum/min while the dendrites on the side arms have a slower 
average velocity 11-12 /im/min even though this is along the 
fast direction. If we slice our data more finely we can see that 
many of the dendrites in fact show an accelerating growth rate. 
This increasing growth rate as the dendrites become more pro¬ 
nounced suggests that the growth is dependent on the local 
available copper surface, which we discuss in more detail in 
the next section. A similar analysis was performed for several 
different samples. All samples were produced using a similar 
low pressure CVD growth, but with small variations on the 
conditions. From the Raman maps of these samples we can 
extract the area and perimeter of the graphlocon as a function 
of time. 

We consider the scaling exponent a defined by A ~ P^, 
where A is the area and P the perimeter. A compact hexagon 
would give a = 2, whereas for a more dendritic shape the 
value should be lower |[35l . In previous works, various scaling 
exponents have been extracted, ranging between 1.43 and 2 
for different flakes under different growth conditions and par¬ 
ticularly sensitive to growth temperature (271. In these works, 
a could only be extracted once at the end of the growth for 
every grown crystal. Here, on the other hand, we are able 
to determine how a evolves during the growth process. Re¬ 
markably, we see that the scaling exponent doesnt change as a 
function of time as shown below. This means that the growth 
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FIG. 3. a) Effect of methane pressure on growth speed, b) Magnitude of the growth velocity shows an increased rate of growth along the six 
primary arms of the graphene crystal, c) Colour map indicating the spectral intensity as a function of Raman shift and position, taken along 
the fast and slow profiles of the crystal as indicated in b. 
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FIG. 4. a) Extracting the average growth rate of the primary arm and several dendrites of a large graphlocon b) Perimeter and area of the 
graphlocon as a function of time for several samples prepared with different growth conditions. Data points are offset on the y-axis to show 
the difference in slope, such that the area of samples D, E and B are scaled by factors of 8,4 and 2 respectively. 


conditions determines a, which is time-invariant, while the 
growth time simply determines the final size of the crystal. 
Hence, for a dendritic crystal, the growth rate in the fastest 
direction is not constant, but a is. 


In figure we show the area versus perimeter of 4 dif¬ 
ferent graphene crystals. A powerlaw fit gives values of a 
ranging from 1.4 to 1.98 for the scaling exponent. The range 
of validity of a single powerlaw extends over the entire size 
of the crystal, demonstrating the time invariance of a during 
the growth process. We also note that the fractal dimension 
increases with an increasing ratio of hydrogen to methane dur¬ 
ing the growth. In one case (sample D) we obtain a value of a 
considerably larger than for the other samples. This sample, 
also shown in figure [^, was prepared with a much larger hy¬ 
drogen flow of 80 seem when compared to all the other sam¬ 
ples, which ranged from 4 to 12 seem H 2 and leads to fewer 
dendrites and a 2. 


Diffusion growth area 

Earlier studies have already examined the time dependence 
of dendrite growth and merging in graphene and found that 
growth is faster at the tip of the dendrite and slower for regions 
in between dendrites and slower at the end of the growth l^ . 
However these studies lack the time resolution to conclusively 
determine the factors affecting the growth rate. To shed light 
on this process, we show that the growth rate correlates with 
the locally available copper surface and argue that the den¬ 
dritic growth is limited by the diffusion of the adsorbed C 
species on the available copper surface. We estimate the re¬ 
lated diffusion rate to be 6 /im as discussed below. 

To gain insight on the locally available copper surface and 
the growth rate we extracted the growth rate along many dif¬ 
ferent profiles along the edge by focusing on the fastest growth 
direction. From the high resolution raman map we can ob¬ 
tain for any given time t, the exact shape of the graphlocon 
as well as the growth rate in any direction. In figure we 
see two examples illustrating results for different dendrites. 
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FIG. 5. Top: Example of an accelerating dendrite. Filled green circle indicates the free copper surface , Sf, while the arrow indicates the 
magnitude of the growth velocity, v. Middle: Growth velocity as a function of time, with large circles indicating the snapshots shown. Bottom: 
an example of a decelerating dendrite. A full movie can be seen in the supplementary material 


The top line shows the behaviour of a dendrite growing into 
free space. As the dendrite grows the shape becomes more 
pronounced and the available copper surface increases. This 
results in an increasing growth rate. The second line shows a 
dendrite growing into occupied space. In this case as the den¬ 
drite grows the available copper surface is decreasing and the 
growth rate is limited by the lack of available surface. This 
results in a decreasing growth rate. The correlation between 
the locally available copper surface and the growth rate was 
averaged across 20 dendrites. We observe an approximately 
linear relationship above 30% available surface and a roughly 
constant growth rate of 5/im/min below 30%. 

In order to determine the approximate radius of the diffu¬ 
sion area we first assume that the growth rate is a linear func¬ 
tion of free surface area over the range of 30%-70%. We fit 
the data in that range considering radii from 1 to 20 /im and 
extract the coefficient of determination for the linear re¬ 
gression. The radius corresponding to the maximum de¬ 
termines our diffusion radius, which we find to be approxi¬ 
mately 6 /im. This value is similar to the observed spacing of 
dendrites seen in figure |^. The mechanism emerging from 
this, is that the methane is adsorbed on the copper surface and 


then diffuses on that surface. Once a graphene nucleation has 
started, the edge of the crystal requires more methane while 
at the same time expelling the accumulated hydrogen from 
the copper induced catalytic conversion of methane to carbon. 
The rate of this reaction depends on the local concentration of 
hydrogen and methane gases. The more methane is available 
within its diffusion area, the faster the growth of graphene. 


CONCLUSIONS 

To conclude, we have developed a technique, which allows 
to image dynamically the growth of graphene mono-crystals 
in addition for providing a tool to label different parts of the 
crystal. This imaging technique enables the assessment of the 
various anisotropic and non-monotonous growth rates respon¬ 
sible for the dendritic and fractal growth of graphene. These 
graphene monocrystals grow from nucleation sites, which are 
typically impurities or deformations in the copper substrate. 
We observed that the graphene nucleation is time dependent 
with individual crystals nucleating at different times indepen¬ 
dent of other factors. Smaller single layer graphene crystals 
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form as a result of delayed nucleation, but the mechanism and 
growth rate are the same for both smaller and larger crystals. 
Finally, the dynamic imaging allows us to draw a precise pic¬ 
ture of the CVD growth process, including the relevant diffu¬ 
sion area of the adsorbed species responsible for the growth. 
Moreover, this technique can easily be expanded to other ma¬ 
terials whenever source elements with different isotopes are 
available, including common elements such as silicon and ni¬ 
trogen. 
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